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Abstract 

Coherently oscillating scalar condensates formed along flat directions of the 
MSSM scalar potential are unstable with respect to spatial perturbations if the 
potential is flatter than c/)^, resulting in the formation of non-topological solitons 
such as Q-balls. Using the renormalization group we calculate the corrections 
to the 0^ potential for a range of flat directions and show that unstable conden- 
sates are a generic feature of the MSSM. Exceptions arise for an experimentally 
testable range of stop and gluino masses when there are large admixtures of 
stops in the flat direction scalar. 
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1 Introduction 



The scalar potential of the Minimal Supersymmetric Standard Model (MSSM) is 
a complicated function of 45 complex squark and slepton fields and 4 complex Higgs 
fields. A natural feature of this potential is the existence of fiat directions, corre- 
sponding to linear combinations of these fields such that there are no renormalizable 
contributions to the scalar potential along the fiat direction beyond the soft SUSY 
breaking terms [Q]. These fiat directions may play a fundamental role in the cosmol- 
ogy of the MSSM, in particular as a natural source of the baryon asymmetry of the 
Universe via the Affleck-Dine mechanism ||^, in which a baryon asymmetry is induced 
in a coherently oscillating condensate of squarks and sleptons. It has recently been 
realized that the cosmology of fiat directions of the MSSM and Affleck-Dine baryo- 
genesis may be more complicated than previously thought |^, |^, |^, ||. If the fiat 
direction scalar potential increases less rapidly than 0^, where is the scalar field 
along the fiat direction ("Affleck-Dine (AD) scalar"), then once coherent oscillations 
of the field begin, the condensate will have a negative pressure making it unstable with 
respect to spatial perturbations 0. These spatial perturbations, which generally 
arise during infiation as a result of quantum fiuctuations of the AD scalar § |^, 
will grow and go non-linear, fragmenting the condensate into lumps which eventually 
evolve into Q-balls of baryon number ("B-balls") [§|, 0) 3- The formation of Q-balls 
has recently been verified by lattice simulations ||10|| . 

The cosmological evolution of the resulting Q-balls depends on the form of SUSY 
breaking. Q-balls from the AD condensate were first discussed in the context of gauge- 
mediated SUSY breaking p, ^ , in which the fiatness of the potential occurs above the 
mass of the messenger fields. It was later realized that Q-balls can also form in the more 
conventional gravity- mediated SUSY breaking models []^ |^. In this case the potential 
is likely to be fiatter than 0^ as a result of radiative corrections. The resulting Q- 
balls are not stable but can be very long-lived, decaying after the electroweak phase 
transition has occured |]^ . As a result Q-balls can protect the baryon asymmetry from 
the effect of lepton number violating interactions combined with sphaleron processes 
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and, if they decay at a low enough temperature, can also act as a common source 
of baryons and dark matter neutralinos, so explaining the similarity of the number 
densities of baryons and dark matter particles when the dark matter particles have 



masses of the order of m^y 11]. They can also enhance the isocurvature density 
perturbations expected from AD baryogenesis in particular inflation models [p!2| , Q . 

The existence and cosmology of Q-balls in the MSSM with gravity-mediated SUSY 
breaking is dependent upon the details of the radiative correction to the flat direction 
scalar potential. The correction must be negative in sign relative to the SUSY breaking 
mass squared term in order for Q-balls to form, whilst the binding energy of the charges 
within the Q-ball and so the total charge and lifetime of the Q-ball depends upon the 
magnitude of the correction Q , as does the fraction of the total baryon number initially 
trapped within the Q-balls In this letter we consider in detail the radiatively 
corrected scalar potential for a range of flat directions using the renormalization group. 

We will show that condensate collapse and Q-ball formation are almost a general 
feature of all the MSSM flat directions, with exceptions only for the cases of d=4 H^L- 
direction and directions with large admixtures of stop. The existence of instabilities 
along the latter directions depend on the mass of the stop in a testable way, as will be 
discussed in the following. 

F- and D-flat directions of the MSSM have been classified and listed in P]. For 
gravity-mediated SUSY breaking the scalar potential along a flat direction has the 
form 0, H 

where m is the conventional gravity-mediated soft SUSY breaking scalar mass term 
(m ^ 100 GeV), is a parameter which depends on the flat direction, and d is the 
dimension of the non-renormalizable term in the superpotential which first lifts the 
degeneracy of the flat direction; it is of the form / = AM^~'^$i ■ ■ ■ (^d- We assume that 
the natural scale of the non-renormalizable terms is M*, where = MpJ v^Svr is the 
supergravity mass scale. 

The equation of state for a field oscillating in a potential 0'^ reads as p = 
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— l)p so that Eq. ([l|) gives rise to the equation of state p = ^p. With < 
the pressure is negative and hence the AD condensate is unstable. The condensate 
fragments and eventually the lumps will evolve dynamically into the state of lowest 
energy, the Q-ball. If > the condensate is stable and no Q-balls will form. 
Although a negative K should be a generic feature of the MSSM, not all flat directions 
will have negative K for all the values of the MSSM parameters. Moreover, as the 
actual value of K dictates the dynamical evolution of the AD condensate and its 
fragmentation, it is of great interest to find out the precise value of K for a given fiat 
direction and for a range of the MSSM parameter values. 

K can be computed from the RG equations, which to one loop have the form 

^ = E + E ^^(E + , (2) 

where a^g and hij are constants, nig is the gaugino mass, A is the A-term, ha the 
Yukawa coupling, and t = \nMx/p- The full RG equations are listed in and we do 
not reproduce them here. We assume unification at t = and neglect all other Yukawa 
couplings except the top Yukawa ht{Mw) = 1 (for definiteness, we choose tanP = 1). 
We shall use quantities scaled by m and denote rrig/m at t = by The potential 
along the fiat direction is then characterized by the amount of stop mixture (where 
appropriate), the values of C, and A, and in the special case of the d=4 HuL -direction, 
on the HuHd -mixing mass parameter pn- 

The mass of the AD scalar is the sum of the masses of the squark and slepton 
fields (pi constituting the flat direction, m| = J^aPf^l ; where pi is the projection of 
(j) along (pi, and = 1- The parameter K is then given simply by 

dmi 



K 



''S 



dt 



(3) 



To compute K, we have to choose the scale p. The appropriate scale is given by the 
value of the AD fleld when it flrst begins to oscillate at if ~ m. Let the value of the 
mean fleld be 0o; the value of K at this scale then determines whether the condensate is 
unstable or not. We may compute 0o from Eq. (|l|) by ignoring the radiative correction 
(i.e. setting effectively = 0) and minimizing the U(l) symmetric part of U (or 



neglecting the A-term in Eq. (|I])). One then finds 



fl = 100 

where in what follows we assume for simplicity that A = 1. 



Table 1. The fiat directions considered 



direction 


dimension 


mass^ 


HuL 


4 


\{ml + III + ml) 


uude 


4 


\{ml^ + ml^ + +ml)] {i^ j) 


QQQL 


4 




{uddy 


6 




{QLdf 


6 





The fiat directions we have studied are listed in Table 1. There is one purely lep- 
tonic direction, H^L, a purely baryonic one involving only the squarks, and directions 
which have both squark and lepton fields. 

In Fig. 1 we show the contours of K for the d=4 uude and QQQL directions in 
the (^,0 -plane; Fig. 2 is for the d=6 {uddY and {QLdY directions. These should be 
representative of all the other directions, too, except for HuL. For ^ ~ typical 
value for K is found to be about —0.05. For all the squark directions with no stop, as 
long as hb and hu can be neglected, K is always negative, and the contours of equal 
K do not depend on A. This is evident from the RGEs Eq. (|^). However, as far as 
flat directions are concerned, all squarks are equal, and having no stop mixture would 
appear rather unnatural. Therefore we have considered the effect of stop mixing in 
the squark directions, which results in A-dependence of the i^-contours, as is depicted 
in Figs. 1 and 2. In the presence of stop mixing if < is no longer automatic even in 
the purely squark directions. The more there is stop, the larger value of ^ is required 
for K < 0. However, even for pure stop directions, positive K is typically obtained 
only for relatively light gaugino masses with ^ < 0.5. In d=4 directions the effect of 
stop mixture is less pronounced than in the d=6 directions, as can be seen from Figs. 
1 and 2. 
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Figure 1: Contours of K for two d=4 flat directions: (a) X = 0; (b) K = -0.01; (c) 
K = -0.05; (d) K = -0.1. The directions are (i) QaQsQL; (ii) QQQL, no stop; (iii) 
u^ude; (iv) uude with equal weight for all w-squarks. 
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Figure 2: Contours of K for two d=6 fiat directions: (a) K = 0; (b) K ^ -0.01; (c) 
K = —0.05; (d) K = —0.1. Tlie directions are (i) QsLd; (ii) QLd, no stop; (iii) u^dd; 
(iv) udd with equal weight for all ti-squarks. 
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In contrast to the squark directions, K was found to be always positive in the H^L- 
direction. This is due to the fact H^L does not involve strong interactions which in 
other directions are mainly responsible for the decrease of the running scalar masses. 
The value of /ih was chosen in such a way that for each value of A and ^, electro weak 
symmetry breaking is obtained at the scale Mw- 

Except for the ifuL-direction, the instabihty of the AD condensate is thus seen to 
depend on the amount of stop mixture in the fiat direction. Since it would be natural 
to expect roughly equal mixtures of the different squark flavours in the fiat direction 
scalar, this means that in principle instabilities, and hence Q-balls in the case of gravity 
mediated susy breaking, could be ruled out experimentally at LHC by measuring the 
mass of the gluino and the stop. To illustrate this, in Fig. 3 we show the the domains 
of positive and negative K in the region of (m^, rrig) -plane corresponding to the range 
—3 < A < 3 and < ^ < 2.5. The Figure is for the {uddY direction with equal weight 
for all H-squarks. For most part K > Q and K <Q regions can be separated, although 
there is a small area below the upper K — Q contour where both values can be found. 
For a fixed m, the K — contour has endpoints which correspond to A — ±3; if one 
were to allow for a wider range in A, this would spread the region between the dashed 
lines towards the lower right-handed corner. Changing the value of m would redefine 
the physical mass scale by a factor m/100 GeV. Thus measuring mi and rrig would not 
alone be sufficient to determine the existence of instabilities. In addition, one needs 
the values of m and tan(5 which naturally will be measured if supersymmetry will be 
found. Very roughly, instability is found when rrig > mi, although the exact condition 
should be checked case by case. 

In conclusion, we have shown that the MSSM scalar condensates in all but the 
H^L flat direction are unstable for a large part of the parameter space. Therefore the 
existence of Q-balls is a generic feature in all the models that incorporate both the 
MSSM and inflation. Moreover, as the existence of instabilities can in principle be 
ruled out by measuring the mass of the stop and the gluino, one may soon be able to 
subject the Affleck-Dine scenario to a direct test. 
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Figure 3: The regions of positive and negative K in the (m^, mg)-plane in the interval 
— 3 < A < 3 and < ^ < 2.5 (between the dashed hnes) for the zidci-direction with 
equal weight for all li-squarks. The solid hue is the K — Q contour, horizontally 
hatched area is where X < 0, vertically hatched area is where X > 0, and there is a 
small region below the upper K = Q -line where both K > Q and K <Q can be found. 
The units are (m/100 GeV) GeV. 
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